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ABSTRACT 

We present approximate formulae describing the changes in neutron star parameters 
caused by the first-order phase transition to an "exotic" state (pion or kaon condensate, 
quark matter) resulting in formation of a mixed-phase core. The analytical formulae 
for the changes in radius, moment of inertia and the amount of energy released during 
the core-quake are derived using the theory of linear response of stellar structure to 
the core transformation. Numerical coefficients in these formulae are obtained for two 
realistic EOSs of dense matter. The problem of nucleation of the exotic phase as well as 
possible astrophysical scenarios leading to a core-quake phenomenon is also discussed. 
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1 INTRODUCTION 

After the discovery of neutron stars it became immediately 
clear that they are unique playgrounds to study the prop- 
erties of super-dense matter in the most extreme physical 
conditions. On the other hand, the very existence of neu- 
tron stars, and their participation in various astrophysical 
phenomena, represented a major challenge for theorists. In 
order to construct neutron star models, one needs to know 
the equation of state (EOS) and other properties of matter 
at densities > 10 15 g cm -3 . 

One of the most intriguing predictions of some theories 
of dense matter is a possibility of a phase transition into 
an "exotic" state. Several phase transitions are predicted 
by dense-matter theories, including pion and kaon conden- 
sation, and de-confinement of quarks (for review see e.g. 
IWeberlll999h . The most interesting case, from the point of 
view of its consequences for neutron star structure, corre- 
sponds to a phase transition which is of the first-order type. 
In such a case, equilibrium phase transition from the nor- 
mal, lower density phase to the pure exotic one, occurs at a 
well defined pressure, Po, and is accompanied by a density 
jump at the phase interface. 

A first order phase transition allows for a metastability 
of the pure normal phase with pressure P > Po- Therefore, 
a metastable core could form during neutron-star evolution 
(accretion, spin down), and nucleation of the exotic phase 
would then lead to formation of a new-phase core, accom- 
panied by a core-quake, energy release, and other phenom- 
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ena, such as shrinking of the radius and a speed up of ro- 
tation. Usually, the radius of the new-phase core is signifi- 
cantly smaller than the original stellar radius, and the effect 
of core-quake can be described within the linear r e sponse 
the ory developed bv [HaerrsekZdunik fc Schaefferl jl986fl 
and IZdunik. Haensel fc Scha effer] (|l 987fl : an earlier Newto- 
nian t heory was presented bv Ischaeffer. Haensel fc Zdunikl 
(1983). Expressions for the changes of neutron-star radius, 
moment of inertia, and energy release were obtained in these 
papers in terms of expansions in powers of the radius of the 
new-phase core, the numerical coefficients of the expansions 
being determined by the EOS of normal phase, the mass of 
the metastable configuration and the density jump at the 
phase transition. 

Relaxing the condition of local (microscopic) electrical 
neutrality opens a possibility of coexistence of two phases of 
dense matter, within a range of pressures, in the form of a 
mixture of the two phases, lower-density (normal - N) and 
higher- density (superdense - S), each of them charged, and 
the mixture being electr ically neutral on the average only 
iGlendenninall99fl . ll992ll . The volume fraction occupied by 
the higher-density phase grows from zero at the lower pres- 
sure boundary, P^ m \ of the mixed phase up to one at the 
upper pressure boundary, Pg™' 1 . If the surface tension at the 
N-S phase interface is not too large, then the mixed phase 
is preferred over a pure phase state (the latter is: N-phase 
for P < P and S-phase for P > P ). 

In this paper paper we study the changes of 
the neutron-star parameters implied by the appearance 
of an exotic phase which forms a mixed-phase core. 
Some astrophysical scenarios which could lead to for- 
mation of a mixed-pha se core were already considered 
in literature (see e.g. iGlendenning. Pei fc Weberl Il997l : 
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Heiselberg fc Hiorth-Jenseiil Il998t IChubarian et alJ l200Ct 
Sp"^^^gS^reio^a^l200a . The phase transition in stellar 



core results from the increase of the central pressure due to 
the pulsar slow-down or mass accretion on the neutron star 
in a binary system. In all these papers the authors assume 
that matter remains in equilibrium during neutron star evo- 
lution, and therefore they use an EOS of matter in ground 
state. However, in a first order phase transition, the new 
phase can appear only via nucleation. Therefore, the star 
acquires first a metastable core of normal phase, in which an 
exotic phase nucleates. Nucleation destabilizes stellar config- 
urations, implies a core-quake, and finally new equilibrium 
is reached with a core of a mixed or a pure exotic phase, the 
actual outcome depending on the detailed kinetics of the 
first-order phase transition. 

In the present paper we derive the formulae for the 
changes of neutron star parameters resulting from the forma- 
tion of a mixed-phase core. These changes are proportional 
to the specific powers of the ratio of the mixed-phase core 
radius to the radius of the last stable pure N-phase config- 
uration. We neglect the effects of rotation on neutron-star 
structure. 

The plan of the article is as follows. In Sect.|5J our nota- 
tion, construction and description of the EOS with a mixed 
phase segment is presented. General lowest-order formulae 
for the changes in neutron-star parameters are derived in 
Sect.|H] Numerical calculations are performed for two recent 
realistic EOSs of the normal phase, in which the exotic dense 
phase is assumed to nucleate. Results for the linear response 
parameters are presented in Sect. where we also provide 
with some numerical estimates of expected effects. In Sect.|5J 
we briefly describe the astrophysical scenarios which could 
lead to a nucleation of an exotic phase and a core-quake. 
Some problems connected with nucleation and formation of 
a mixed-phase core are discussed in Sect.|S| Finally, in Sect. 
Qwe summarize our results and we also discuss their possi- 
ble applications, and perspectives of their generalization to 
rotating neutron stars. 



2 EOS WITH A MIXED-PHASE REGION 

Let us consider a general case of a first-order phase transition 
between the N and S pha ses of dense matter. As shown by 
iGlendennind dl99ll Il992t) , relaxing the microscopic charge- 
neutrality condition can make a mixed-phase state energet- 
ically preferred provided the surface tension and Coulomb 
contributions are sufficiently small. 

Assume the thermodynamic equilibrium of a multi- 
component and multi-phase dense matter, neglecting the 
Coulomb and surface contributions to the thermodynamic 
quantities. The elementary constituents of the matter are 
hadrons (h), which may be baryons, quarks, and strongly- 
interacting meson (pion and kaon) condensates, as well as 
leptons (electrons and muons). The energy densities in both 
phases depend on the number densities of the matter con- 
stituents in these phases, 



£ = £ ({n h },n e ,n M ) , £ =£ ({n fe }, n e , n„ 



(1) 



As the translational invariance may be broken within 
a phase, the number densities are actually the volume- 
averaged ones. We assume that the size of the region occu- 



pied by a non-uniform phase is sufficiently large compared to 
the characteristic length-scale of the non-uniformity. Then 
the volume averages within each phase are well defined. 

The corresponding electric-charge densities (in the units 
of the elementary charge) and baryon densities (in the units 
of nucleon baryon charge) are given by 
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where qn and bh are, respectively, the electric and baryonic 
charges of a hadron h. 

The thermodynamic equilibrium of a mixture of phases 
N and S, at a fixed average baryon density rib, will be cal- 
culated by minimizing the average energy density 



£=(1~ X )£ N + X £ S ■ 
under the condition 

/i \ N . S 

n b = (l- X)n h + xn-b 



(3) 



(4) 



and under the constraint of average (macroscopic) electrical 
neutrality 



Pd. = (1 - x)Pei. + XPel. = 0. 



(5) 



For example, let us consider a first-order phase tran- 
sition implied by kaon condensation. Let us assume that 
the characteristic length-scale of the hadron electric-charge 
inhomogeneities is much smaller than the electron and 
muon screening length s (this is not always true, see e.g. 
iNorsen fc Rcddv 2001J). Then the electron and muon den- 
sities can be considered as uniform, 

N S _ N S i a \ 

n e = n e = n e , ra M = n M = n M . (6) 

We have to determine the values of eight variables (four nu- 
cleon densities, two lepton densities, kaon density, and the 
volume fraction x) by minimizing £ , Eq. 0, under the con- 
ditions given by Eqs. J1J and JSJ. This leads to a set of non- 
linear equations relating the thermodynamic variables, each 
of these equations having a clear physical meaning. Mechan- 
ical equilibrium between the two phases requires (we remind 
that surface and Coulomb contributions are neglected) 
pN = p s ^ 

The strong interactions imply the equality of chemical po- 
tentials of nucleons in the two phases, 

,,N S _ N_S_ /on 

Pn — Pn — P™ 1 pp — pp — Pp> \°J 

Finally, weak interactions involving hadrons and leptons 
lead to 



pn — Pp T" Pe , Pe — pK~ 



Pyi = Pe 



(9) 



Together with Eqs. @ and JSJ, we get eight equations for 
eight thermodynamic variables. The solution corresponds to 
the thermodynamic equilibrium at a fixed rib and under the 
constraint of macroscopic electrical neutrality. Let us first 
discuss the character of the bulk equilibrium as a function 
of the pressure, Fig. For P < ■ , the equilibrium 



is realized by the pure N-phase. For P T 



(m) 



< P < P 



(m) 



the equilibrium corresponds to a mixed m-state of both 
phases (the m-phase). The volume fraction occupied by 
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Figure 1. Baryon chemical potential /j b versus pressure P in the 
presence of an equilibrium first order phase transition, at P = Pn, 
between the N and S phases. The N phase, if stable with respect to 
the phase transition into the S-phase, is represented by the solid 
line; its dotted continuation corresponds to an over-compressed 
state, mctastable with respect to the transformation into the S 
phase. Analogous notation is used for the S phase, where the dot- 
ted continuation of the S-line corresponds to a metastable under- 



compressed state. The mixed 
represented by the dash-dotted line 
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Figure 2. Same as in Fig. HI but for P versus The mixed 
phase (pX < P < Pq ) is represented by the dash-dotted line. 
For further explanations see the text. 



the S-phase increases monotonously with P, from zero at 



P = P 



(m) 



to one at P — P c 



(m) 



For P > Pg we have 



a pure S-phase. The calculations of a mixed-phase state, 
where kaon-condensed m atter coexisted with baryon mat- 
ter, w e re pe rforme d bv IGlend enning fc Scha ffner-Bielichl 
<199St 119991) and iNorsen fc Reddvl feOOll) . The mo" 
els of dense matter with a mixed phase of de-confined 
quark matter coexisting with baryon pha se were con 

and 



structed bv |Heiselbere . Pethick fc Staubol 1^93) and 
IGlendenning fc Peil <ll995f) . The importance of the phase- 
interface (surface tension and curvature energy) effects 



for forming the mixed-pha se state and for its spatial 
struct ure was studied bv iHe iselberg. Pethick fc Staubol 
Jl995t: IChristiansen fc Glendennind <1997tl : 

Christiansen, Glendcnning fc Schaffner-Bielichl (2000) 
and lNoreai^^Rfid^n200lTr ~A mixed-phase state affects 
the EOS of dense matter as visualized in Fig. For the 
sake of comparison, we show also a standard first-order 
transition between the pure N and S phases. For a pure 
N-S phase transition, the densities «n < rib < ns could not 
exist in the stellar interior because P should be monotonous 
there to create an outward directed force which balances 
at each point the gravitational pull. On the contrary, the 



mixed-phase layer of density n 



(m) 



< 7ib < n. 



(m) 



can well 



exist in the stellar interior, with the pressure increasing 



from PjJ, m ' at the top of the layer to Pg m; at its bottom. 
Mixing of the N and S phases softens the EOS as compared 
to the pure N-phase EOS, but the effective softening is 
weaker than in the limiting case of a pure-phase transition 
between the N and S phases. 

The surface and Coulomb effects bring positive con- 
tributions to £ and Hb- They affect the size and shape of 
the structures within the mixed-phase layer. For a periodic 
struc ture, the virial theorem (see, e.g.. |Pethick fc Ravenhaiil 
Il995l and references therein) tells us that the surface contri- 
bution is twice the Coulomb one. Both are pushing up the 
value of At b m '(P)- They are particularly important at the 
edges of the mixed-phase region, where the droplets of one 
phase within the dominating one are small. It is clear that 



(m) 



these effects increase P^™' 



and decrease Pg™', narrowing 



the mixed-phase layer in the stellar interior. If sufficiently 
large, the surface and Coulomb contributions can entirely 
remove the mixed phase: the difference in fib(P) of the pure 
and mixed phases is usually small. It will occur if surface 
tension is greater than some critical value (a > a" cr it) so 
that ;i m) fP ) > t£(P) for P < P n and ^. m) (P) > Mb(P) 



for P > Pn. IHeiselberg. Pethick fc Staubol J1995T) obtained 
(Tcrit — 70 MeV fm~ 2 for a transition from nucleon (N) to 
quark (S) matter. The actual value of the surface tension 
for the quark-matter droplets in baryonic medium is very 
poorly known, a = (10 - 100) MeV fm" 2 . 



3 CHANGES OF THE STELLAR 

PARAMETERS FOR THE MIXED-PHASE 
TRANSITION - LINEAR APPROACH 

We assume that at a central pressure P c = P cr it the nucle- 
ation of the S-phase in a super- compressed core, of radius 7"n> 
of configuration C, initiates the phase transition and forma- 
tion of mixed-phase core of radius r m in a new configuration 
C* , as presented on Fig. El Transition to a mixed phase oc- 
curring at r m is associated with a substantial drop in the adi- 
abatic index of matter, defined as 7 = (nb/P)dP/drib, from 
7n to 7 m : mixed phase is much softer than the pure one. This 
is because in the mixed phase the increase of mean density 
is reached partly via conversion of a less dense N phase into 
denser S phase, and therefore requires less pressure than for 
a pure phase. The effect exists for any fraction of the S phase, 
even in the limit of the fraction of the S-phase tending to 
zero, and leads to a discontinuity of 7 at p^ 11 ' = p m . An ex- 
ample of a very dramatic drop of 7 at p m is giv en in Fig. 20 of 
lAkmal. Pandharipande fc Ravenhh all (199q). Other exam- 
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Figure 3. Central pressure P c versus central matter density p c for 
neutron-star configurations based on a pure N-phase EOS and an 
EOS with a mixed-phase segment. Solid line denotes stable states, 
dash-dotted line - the states which a meta-stable with respect 
to the transition to a mixed-phase state. For a critical central 
density p C rit the S-phase nucleates in the super-compressed core 

of configuration C, this results in a transition C ► C* into a 

stable configuration with a mixed-phase core and central density 
p*. Both configurations C and C* have the same baryon number 
A. 

pies are given bv lPons et alJ i2000L see Fig. 1 in their article), 
where the transition from the solid line (no kaon condensate) 
to dashed line representing mixed phase is clearly accompa- 
nied by a drop in the derivative of pressure with respect to 
density. In Fig. 0] we show two exa mples of realistic EOSs 
with a mixed phase segm ent, one bv |Ofende nningl l)l997i Ta- 
ble 9.1) and the other bv lPons et alJ (1200(1 GM-GS-U120'). 
As one can see, near the mixed-phase transition point the 
behaviour of the logarithm of pressure is linear with respect 
to the logarithm of baryon density; it clearly means that the 
polytropic approximation is valid in the small mixed-phase 
core regime. This is precisely the case we are interested in - 
we will henceforth benefit from expansion in powers of the 
core radius r m . The validity of the 7 m = const, approxima- 
tion in the limit of small r m can be also phrased differently, 
using the s mall-r m expansions. Using techniqu es described 
in detail in lZdunik. Haensel fc Schaefferl il987t) . we can see 
that the difference between the value of 7 at p m , and at the 
center of the C* configuration in Fig. is quadratic in r m , 
so that 7* = 7 (p*) — 7m + 0(r,n)- As we show below, the 
leading changes due to a core-quake are proportional to the 
fifth (radius, moment of inertia) and seventh (gravitational 
mass) power of r m . They are obtained by putting 7* ~ 7 m . 
Inclusion of quadratic correction in the expression for 7* 
contributes to the higher order terms, i.e., seventh power of 
r m in the case of radius and moment of inertia, and ninth 
power of r m in the case of gravitational mass. 

Having a pair of EOSs, one with and one without a 
mixed-phase softening, we compare the hydrostatic equilib- 
ria of neutron stars corresponding to each of these EOSs. 
The models which we calculate are non-rotating, spher- 
ically symmetric solutions of Einstein's equations, which 
in this special cas e reduce to the Tolman-Oppenheimer- 
Volkoff equations jTolmanl Il939t lOppenheimer fc Volkofll 
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Figure 4. Upper panel: logarithm of pressure P versus the loga- 
rithm of baryon density f or sample mixed-ph ase EOS (quark 
matter - baryon matter) from lGle ndcnning ( 1997), Table 9.1. The 
filled circles correspond to the tabulated points. The discontinu- 
ous drop in 7 at the N-m interface is clearly seen. Notice that the 
deviation from the linearity of the EOS of the shown m-phase seg- 
ment in the log(njj)— logP plane, is small. This means, that within 
the considered density interval, which extends from 0.26 fm — 3 to 
0.35 fm — 3 , polytropic approximation for the m-phase is good. 
Lower panel: similar plot for the GM-GS model with kaon opti- 
cal potential = —120 MeV, producing the EOS with mixed 
phase segment (k aon-condensate - baryon matter), obtained by 
iPons et alJ l2000l) . The filled circles correspond to the tabulated 
points. The polytropic approximation for the mixed-phase with 
kaon condensate appears to be very good. 

Il939l) . The neutron-star models are labeled by the central 
density p c . The configurations calculated with these two 
EOSs are identical up to p c = Pn°'; the configuration with 
such central density is denoted by Co, and will be treated as 
a "reference configuration". 

When the central density exceeds pjj" , the models 
based on these EOSs are different, due to the appearance 
of a softer mixed-phase core in configurations corresponding 
to the mixed-phase segment of the EOS. For a fixed baryon 
number A, greater than a baryon number Aq of the reference 
configuration Co, we compare the global parameters, such as 
mass-energy, radius, and moment of inertia. Their difference 
corresponds to the changes in these parameters implied by 
the phase transition in the stellar core. 

From now on we will restrict ourselves to a lin- 
ear response of neutron star to the appearance of the 
mixed phase core. The calculation is based on express- 
ing the change in the density profile, due to the pres- 
ence to a small core, as the combination of two in- 
dependent solutions o f the linearly perturbed equations 
of stellar structure I lIggrrggL Zdunik fc Schaefferl Il98fil : 
IZdunik. Haensel fc Schaefferll987l) .' The presence of a denser 
phase in the core changes the boundary condition at the 
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phase transition pressure P^ and allows us to determine 
the numerical coefficients in the expression for the density 
profile change. The leading term in the perturbation of the 
boundary condition at the edge of the new phase results from 
the mass excess due to the lower stiffness, and higher density 
of the new phase as compared to those of the N-phase. 

Let us compare basic expressions obtained when 
the new c ore consists of a pure S -p hase, consid- 
ered bv lHa ensel. Zdunik fc Schaeffeil lll986l) and 
IZdunik. Haensel fc Schaefferl l)l987h . and in the present case 
of a mixed-phase core. For a pure S-phase core the phase 
transition is accompanied by density jump pn — > ps which 
leads to the lowest-order expression for the core-mass excess 
(with respect to the pure N-phase configuration), 

4 3 

Sm COIC = -7r(ps - Pk)vs + 0(r s ), (10) 

where powers (rs) with I > 3 have been neglected. In the 
case of a mixed-phase core, considered in the present pa- 
per, the lowest-order expression vanishes, because there is 
no density jump at the core radius. 

Let us introduce the notation appropriate for descrip- 
tion of the change of the properties of matter at the N- 
phase - m-phase transition point. The transition occurs at 

P = Pm 



„C«0 



N . Density does not change and is equal 
p m = /% J , while the adiabatic index changes from 7n on 
the N-side to a lower value 7 m on the m-side. 

To the lowest order in r m , the mass excess in our case is 
due to the difference in the stiffness (adiabatic index) of the 
matter in N-phase and m-phase. Expanding term by term, 
wc get 



5m c 



4np„ 
45 



(1 + xm)(l + 3x m )(n 2 m - /sSr)r5, + 0(rl)(ll) 



where x m is the ratio of pressure P m , and the energy density 
PmC 2 at the phase transition point, x m = P m /p m c 2 . The 
parameters are defined by: 



k n = 4TvGp m /v N , K m = 4nGp m /v 



2 

m ! 



(12) 



where vn and v m denote the sound velocity on the both sides 
of the core boundary. Let us remind that sound velocity v 
is related to the adiabatic index by v 2 — r yP/(p + P/c 2 ). 

The main parameter of the linear theory - the radius 
of the core of mixed phase r m is connected with the density 
range of metastability which can be achieved by the N phase 
of matter, i. e., the difference between p C rit an d p^ , via 
relation (see lZdunik. Haensel fc Schaefferlll987l) : 



Ap crit 



Pcrit 



(m) 

Pn 



6 



«N(1 + £m)(l + 3x m )rn 



(13) 



Strictly speaking, the lowest order works only for suf- 
ficiently small value of r m . From Eq. 11 U and Eq. 1121 we 
have 



4 2 
<5m corc = — 7rp m (l + x m ) (1 



3x n 



4vrGpi 

x 7N/7n 



1) ■ r m 5 + 0(rl) 



(14) 



The new phase manifests itself in the boundary condi- 
tions by the presence of the prefactor ("fn/fm— 1), which acts 
similarly as the prefactor (ps/pN — 1) in the case of a first 
order phase transition from pure N-phase to pure S-phase. 



The fact that in our case the density stays continuous, re- 
sults in the linear-response effects which are proportional to 
a power of the core radius greater by two than in the case 
of the transition to a pure S-phase. 

On the basis of the above analytical considerations, we 
expect that the relative changes in stellar parameter Q con- 
tain a prefactor (7N/7m — 1) and are proportional to r m in 
the case of Q = R, I and to r m in the case of mass-energy 
E — Ale 2 . Summarizing, we expect the following form of the 
lowest-order expressions: 



SQ 



Q* -Q 
Qo 



-(in/jm - 1) -Pq ■ (r m )\ 



(15) 



where f m = r m /Ro, and I — 5 for radius R and moment of 
inertia I, and I — 7 for the energy E. The coefficients (3q 
are functionals of the reference configuration, /3q = /3q[Co]- 



4 CHANGES IN STELLAR PARAMETERS 
FOR REALISTIC EOSS 

In order to explore how a realistic neutron star will react 
on the appearance of a mixed-phase cor e, to describe the N- 
phase we used two rece nt EOSs: SLv oflDouchin fc Haensell 
il200ll) and FPS EOS of lPandharipande fc RavenhalHl98Sft 
Both EOSs describe in unified way (i.e., starting from a 
single effective nuclear Hamiltonian) both the crust and the 
core of neutron star. They assume that neutron star core is 
composed of neutrons, protons, electrons and muons. 

In the actual calculations, we approximated the mixed- 
phase softening by replacing the mixed-phase segment of the 
EOS starting at n m = n^ 11 ', by a polytrope with 7 m < 7n- 
As we explained in Sect. 3, this is equivalent to neglect- 
ing the r-dependence of 7 within the mixed-phase core, 
and it does not change the leading terms of the small r m - 
expansions. Also, as shown in Figs. 2] this approximation is 
excellent for not too large mixed-phase core. 

We constructed a large set of EOSs with a mixed-phase 
segment, using many values of the phase-transition density 
n m and several values of j m . In this way we were able to 
study the linear response of neutron star to the appearance 
of a small mixed-phase core for a wide choice of the mixed- 
phase parameters. 

Note again that the coefficients /3q depend only on the 
properties of the reference configuration Co - the EOS of 
the mixed-phase (here approximated by a polytrope) enters 
merely via 7 m in the factor (7r-j/7m — 1) in Eq. 1151 . 

We apply the arguments presented in Sect.|3]to the EOS 
constructed in the manner described above. By changing the 
value of pm, we obtain a family of reference configurations 
Co, and by comparing C and C* we extract the values of the 
linear response coefficients /3q[Co]- 



4.1 Results for SLy EOS of the N-phase 

The values of (3q [Co] are plotted in Fig. An important 
parameter relevant for the linear response to a perturbation 
of equilibrium configuration is the adiabatic index of the N- 
phase EOS. We plot its value at the center of the reference 
configuration Co versus M[Co] and its central baryon number 
density in Fig. |S| As it can be seen on Fig. for a wide 
astrophysically interesting range of neutron-star masses, the 
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Figure 5. The linear response parameters /3q, versus the mass 
Mo of the reference configuration Co • The N-phase is described by 
the SLy EOS (solid lines, M max = 2.05 Mq) and, for comparison, 
by FPS EOS (dashed lines, M max = 1.79 M ). 
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Figure 6. The value of the adiabatic index at the center of 
non-rotating configuration for the SLy EOS (upper panel) and 
FPS EOS (lower panel), plotted against the gravitational mass 
(bottom horizontal axis), and baryon number density n c at the 
star center (upper horizontal axis). Vertical dotted line marks the 
maximum allowable mass (M max = 2.05 Mq for the SLy EOS, 
and M max = 1.79 Mq for the FPS EOS). 



values of the functionals /3q[Cq] stay almost constant. The 
plateau values for SLy EOS are 0f at = 2.0, /3^ 1 = 0.8, and 
^piat _ q qq £ Qr momen t of inertia, radius and mass- 
energy, respectively (flj lat = 2.4, /3^ lat = 1.0, and /3g lat = 
0.05 for FPS EOS, see next sub-section for details) 

The values of /3q increase rapidly with decreasing mass 
below 0.4 M . This is due to the speci fic generic fea tures of 
the realistic EOS of neutron star crust l)Haensell200ll ). which 
results in the existence of the minimum neu tron star mass 
Mmi n (see lHaensel. Zdunik fc Douchinll2002l . and references 
therein). For M — > M m i n neutron stars become more and 
more "soft" with respect to the fundamental mode of ra- 
dial perturbation, and become unstable at M m i n . Therefore, 
at the same mass-excess 8m cola the effect of neutron star 
"shrinking" due to an increased gravitational pull by the 
mixed-phase core grows, and very rapidly non-linear effects 
become important. 

The "numerical plateau" values deserve an ad- 

ditional comment. The plateau region extends within 
O.5M <M[C ] <0.9M max . It is easy to verify that our 
plateau values of /3q multiplied by ~ 3 are very close 
to the linear response coefficients for the pure S-phase 
co re, ao, obtained for a medium-stiff EOS of the N-phase 
bv lHaensel. Zdunik fc Schaefferl il986T> . Clearly, apart from 
some normalization constant, the linear response coefficient 
does not depend on the nature of the new-phase core, which 
enters only via (ps/pN — l)(rs) z in the case of a first-order 
phase transition to a pure S-phase, and via (jx/jm — l)(f m ) i 
for a mixed-phase core. This results from the fact that 
the change of the stellar parameters Q depends on the 
mass excess in the core via the relation (leading term): 
8Q ~ (dQ I dm)Sm COIC and the numerical factors multiply- 
ing (7N/7m - l)(Tm) 5 or (ps/pN - lX^s) 3 in 5m corc in the 
case of second and first order transition respectively differ 
by j^k^R 2 which is approximately 1/3. 



4.2 Results for FPS EOS of the N-phase 

The values of /3q [Co] are plotted in Fig. |H| where they can 
be easily compared with those obtained for the SLy EOS. 
The plateau values are /3^ lat — 2.4, ,0^ lat = TO, and /3^ lat = 
0.05. As we see, they depend rather weakly on the EOS 
of the N-phase, provided it is realistic. The Mo range of 
the plateau region, (0.6 — 1.6) Mq, is significantly narrower 
than for the SLy EOS. One can try to point out reasons 
for this. First, on the high-mass side, the difference results 
from the fact that M max = 1.79 Mq, to be compared with 
-A^max — 2.05 M . Otherwise, the quasi-constancy is valid, 
similarly as for the SLy EOS, up to about 0.9 Af max , which 
seems to be generic. "Softness" of hydrostatic equilibrium 
at M close to M max with respect to radial perturbations 
is the general-relativistic feature, and does not depend on 
the EOS of the N-phase. On the contrary, on the low-mass 
side, the difference in the /3q behavior reflects the differences 
between the EOSs at sub-nuclear densities. This difference 
between the SLy a nd FPS EOSs was studied, albeit in a 
different context, bv lHaensel. Zdunik fc Douchinl i2002l) . In 
the relevant stellar mass range, the FPS stars are less bound 
and therefore, with decreasing mass, its /3 coefficients start 
to increase earlier than for the SLy EOS. 
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4.3 Estimating changes in stellar parameters in a 
core-quake 

Let us now calculate the changes in stellar parameters. First, 
the estimates will be obtained for the SLy EOS of the N- 
phase. Consider Co configuration of mass 1.4 M calculated 
for the SLy EOS. It has R = Ro = 11.73 km, I = I = 1.37 x 
10 45 g cm 2 , and a central adiabatic index 7n — 2.94. Using 
the plateau values of the /3q coefficients, we can rewrite the 
formulae for the changes in neutron-star parameters in a 
form suitable for making numerical estimates. For the stellar 
radius we get, 

AT? ~ -4.23 x 1CT 5 • (2.94/ 7m - l)(r m /l km) 5 km. (16) 

Assume 7 m = 1.5. Then, a 1 km core implies shrinking by 
about 4 cm, a very small star-quake indeed, but still larger 
by more than an order of magnitude than that associated 
with macro-glitches in the pulsar timing. However, the rise 
of the star-quake amplitude with r m is very steep and for 
r m = 4 km we get AR — 42 m, rather impressive even by 
the terrestrial standards. 

Expression for the fractional change in I is 

AI/Io ~ -9 x 10" 6 • (2.94/ 7m - l)(r m /l km) 5 . (17) 

Assume as before that 7 m = 1.5. Formation of a 1 km core 
implies speed-up of the neutron star rotation by Afl/fi ~ 
—AI/I ~ 10 -5 , one order of magnitude larger than in the 
biggest pulsar macro-glitches. If the core radius is 4 km, 
then Afi/Q ~ —AI/I ~ 9 x 10~ 3 , a very distinct feature of 
a neutron-star core-quake. 

Let consider now the energy release. We have 

AE ~ 4.1 x 10 45 ■ (2.94/ 7m - l)(r m /l km) 7 erg, (18) 

which means that formation of a 1 km core with 7 m = 1.5 
will release a total energy of 4 x 10 45 erg. Corresponding 
amount of released energy for a core with r m = 2 km and 
r m — 4 km is 5 x 10 47 erg and 6.4 x 10 49 erg, respectively. It 
may be worthwhile to compare it with energies associated 
with a few known astrophysical processes involving neutron 
stars. For example, the energy associated with Soft Gamma 
Repeaters (SCR) outbursts, which are now believed to be 
caused by a extremely strong mag netic field, is ~ 10 44 -10 45 
erg. The Vela pulsar macro-glitches release from 10 38 to 10 42 
erg. The thermonuclear Type I X-ray bursts involve energies 
10 36 - 10 39 erg, but in the case of the so-called super-bursts 
the released energy can be as high as 10 43 . By adjusting the 
core radius, the energies associated with core-quakes can be 
made similar to those released in phenomena listed above. 
It should be however noted that the core radius is not a free 
parameter in our approach but is determined by the range 
of metastability possible in the considered physical situation 
which is rather poorly known and model dependent. The 
leading terms in Eqs. (1161181 for the radius r m = 1 km 
correspond to the degree of metastability, which is usually 
called "over-compression" of the order of 0.5%, i.e Ap cr it ~ 
0.005, see Eq. O . 

Up to now, we performed all numerical estimates for 
the SLy EOS of the N-phase. For the same values of r m 
and 7 m , the values of AT? etc., for the FPS model of the 
N-phase, will be somewhat different, but these differences 
are of no practical importance. What is really crucial, is the 



strong dependence of ATi, AT, and especially of AE, on the 
mixed-phase core radius. 



5 ASTROPHYSICAL SCENARIOS 

One can consider several astrophysical scenarios of the for- 
mation of a mixed-phase core in neutron star. In all cases 
we are dealing with a two step process. The first step con- 
sists in the nucleation of the S-phase in the me tastable 
N-phase medium , stud i ed bv [ H aensel fc Schaeffej il982l) : 
iMuto fc Tatsumil il992ft : llida fc Satol (Il99lll998ft . The sec- 
ond step consists in the relaxation of the system to a stable 
mixed state in which N-phase coexists with a S-phase. 

The first scenario is connected with neutron star birth, 
in which proton-neutron star is formed in gravitational col- 
lapse of a massive stellar core. Central core of a proto- 
neutron star has a relatively low density, compared to central 
density of neutron star in which it will transform eventu- 
ally (after a few minutes). This is due to large lepton frac- 
tion (30-40%) resulting from neutrino trapping. The diffu- 
sion of neutrinos and resulting deleptonization of matter is 
connected with softening of the neutron-star core EOS and 
matter compression to higher density. Simultaneously, the 
core is heated (up to 30 — 60 MeV/fce) by the deposition 
of the energy of strongly degenerate electron neutrinos via 
their down-scattering on the constituents of dense medium. 
The optimal conditions for nucleation of the S-phase are 
just after deleptonization (i.e., ~ 10-20 seconds after the 
bounce ending the collapse phase). The thermal fluctuations 
at T ~ 40 Me V/ fee are expected to be sufficient to overcome 
energy barriers, resulting from the Coulomb and surface ef- 
fects, and to form a mixed-phase core in thermodynamic 
equilibrium. 

Central compression of matter during slow-down of pul- 
sar rotation is a second scenario to be considered. Let the 
pulsar rotation frequency be v. The effect of rotation on 
the pressure distribution is to lowest order cx v 2 . There- 
fore, the rate of the quasi-static increase of central pressure 
is P c oc —vv. The central-core temperature is < 10 9 K, so 
that nucleation can proceed only via quantum fluctuation in 
super-compressed core. Nucleation induces pressure deficit, 
collapse of the central core accompanied by matter heating, 
and energy release. During these violent processes involving 
matter flow, N and S-phase mix forming a mixed-phase core. 

A third scenario involves neutron star accreting matter 
in a close binary system. Central pressure increases due to 
the gravity of accumulated layers of accreted matter, at a 
rate P c oc M. The core temperature remains below 10 9 K 
(see e.g. iMiralda-Escude. Haensel fc Paczvriskilll990l) . and 
therefore quantum fluctuations are crucial for the nucle- 
ation process to start. Then the pressure deficit triggers core 
collapse, accompanied by matter heating and flow, and a 
mixed-phase core forms. 



6 SOME PROBLEMS CONNECTED WITH 
NUCLEATION OF A MIXED 
NORMAL-EXOTIC PHASE 

Let us assume that the phase transition from a pure N-phase 
to a pure S-phase is the first-order one, at P = Po (Fig. 
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0. Let us further assume, that within the pressure interval 
< P < -Pg™' the thermodynamic equilibrium of dense 
matter is realized in the form of the mixed NS phase. 

Consider a neutron star, where the central pressure 
P c increases in time due to a mass accretion in a binary 
system or due to a spin-down of neutron-star rotation. 
The central pressure increases on a characteristic timescale 
icvoi = P c /P c - Let us assume that initially the stellar core 
consisted of a pure N phase. In this case the internal tem- 
perature does not exceed 10 9 K, the thermal fluctuations 
are negligible and the S phase can nucleate only via quan- 
tum fluctuations on a certain timescale r. However, as long 
as P c < Po, the nucleation of the S-phase in the quantum 
regime is impossible (r = oo). The actual nucleation will 
start at some P c = Pnuci > -Po, by a formation of a single 
droplet. If the expansion rate of the first drop is larger than 
the formation rate of other droplets in the core, then the 
pure S phase will fill the central core with P < Po. The S- 
phase core will then coexist with the outer layer of the N 
phase, with a baryon density drop ns — nj$ at the interface. 
Such a scen ario seems likely in the case of the quark-matter 
nucleation dlida fc Satdll997t Il99gfi 

The central temperature of a proto-neutron star stays 
high, T c > 10 11 K, for some ten seconds. This may be 
enough to complete the nucleation of the S phase below 
Po and mix it with the N phase, achieving the thermody- 
namic equilibrium. It is not excluded, however, that the core 
will have no time to evolve to thermodynamic equilibrium 
in the low-temperature limit, because of a too rapid cooling 
to T « 10 10 K. It is therefore possible that the final mixed 
state with some spatial structure will be quite different from 
the strict ground state of the core at T = 0; the mixed phase 
can remain "frozen" in some metastable state in the core or 
its part. 

Finally, let us mention two difficulties in for ming a 
mixe d phase in kaon-condensed neutron-star cores l|Norsenl 
2002). It is difficult to nucleate kaon condensate because of 
the slowness of weak interaction processes; very high tem- 
perature T > 10 11 K and low kaon effective masses are 
required for the condensation due to thermal fluctuations. 
This may happen only in massive newborn neutron stars, 
with particularly high central temperatures and densities, 
where a mixed-phase core could be formed. Medium-mass 
neutron stars have insufficient central density to nucleate 
the kaon condensate at their birth. On the other hand, high- 
mass neutron stars, which gain their mass by accretion, can 
remain in a metastable state forever, because their internal 
temperature is too low for nucleating kaon condensate in 
their cores. 



7 CONCLUSIONS 

Formation of a mixed-phase core implies shrinking of neu- 
tron star radius, spin-up, and is associated with energy re- 
lease. Changes in neutron star parameters strongly depend 
on the size of the mixed-phase core. The decrease of stellar 
radius, AR, and of the moment of inertia, AI, are propor- 
tional to the fifth power of the core radius. The total energy 
released during the core-quake, AE, is proportional to the 
seventh power of the core radius. These powers are higher 
by two than the powers of the core radius in the formulae 



expressing the changes of the stellar parameters implied by 
the formation of a pure high-density phase core. 

Apart of the core radius, the formulae involve also 
the drop of the adiabatic index at the mixed-phase 
core edge, and a numerical coefficient determined by 
the normal-phase configuration. In the stellar mass range 
0.5M o <M<0.9M max the numerical coefficients for AR, 
AI and AE depend rather weakly on M and can be re- 
placed by their "plateau value" . 

This paper completes the study, within the linear 
response theory, of the effect of phase transitions in 
the neutron star core on the stellar structure. It ex- 
tends the methods initiated in the 1980s for a s ingle 
(pure ) dense phase core jS chaeffer. Haensel fc Zdunikl 
1 19831: lHaensel. Zdunik fc Schaefferl Il986t 

IZdTmik" Haensel^^5chaefferT " ^)87l) to the case of a 
mixed-phase core. Our results are presented in the form of 
formulae which are convenient to make numerical estimates, 
and which become very precise for small mixed-phase cores, 
when the brute-force calculations are impossible because of 
limited numerical precision. 

The final state reached in a phase transition in neutron- 
star core depends on the kinetics of this process. Even if 
the mixed state is the stable final configuration, it does not 
mean that it is reached, because a one-phase core may form 
instead. Mixed-state is particularly difficult to be formed 
if a significant strangeness production is necessary in the 
dense-phase nucleation process. At temperatures and den- 
sities characteristic of the cores of accreting neutron stars, 
the mixed-phase core may never form. 
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